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Abstract 

We present the results of a full QCD lattice calculation of the flavor singlet axial 
\ vector coupling G\ of the proton. The simulation has been carried out on a 16 3 x 

32 lattice at j3 = 5.6 with nf = 2 dynamical Wilson fermions. It turns out that 
the statistical quality of the connected contribution to G l A is excellent, whereas the 
disconnected part is accessible but suffers from large statistical fluctuations. Using a 
0^ ' 1st order tadpole improved renormalization constant Za, we estimate G\ = 0.20(12). 



p; 1 Introduction 

The flavor singlet axial vector coupling G\ of the proton, 



s v G a = (P\ui^ 5 u + dryffsd + s-f^ 5 s\P) , (1) 

where s M denotes the components of the proton polarization vector, has been the target 
of an intensive research activity both in experimental and theoretical elementary particle 
physics over recent years. About one decade ago, the European Muon Collaboration (EMC) 
experiment found an unexpectedly small value, G\ = 0.12(17), from the measurement of 
the first moment of the spin dependent proton structure function gl in deep inelastic 
polarized muon-proton scattering |l[]. Since G\ can be interpreted in the naive parton 
model as the fraction of the proton spin carried by the quarks, the EMC result became 
known as the so called 'proton spin crisis'. 

A number of succeeding experiments have been performed in the meantime 0, ^ [|, || [| . 
The latest analysis ||, including proton, neutron, and deuteron data finds 

G\ = 0.29(6) , (2) 

at a renormalization scale /i 2 = 5GeV 2 . This is still far away from the Ellis- Jaffe sum 
rule[0 expectation G\ ~ G\ = 0.579(25) f8j, which is found in the Okubo-Zweig-Iizuka 
(OZI) limit of QCD. 
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Figure 1: Connected (a) and disconnected (b) contributions to the axial vector density 
amplitude of a nucleon. Please note that gluon lines connecting the quark lines are omitted. 

As has been pointed out by Veneziano H a possible nontrivial topological structure 
of the QCD vacuum, which induces a nonzero contribution to the divergence of the axial 
vector flavor singlet current via the axial anomaly, might induce a sizeable effect on the 
value of the matrix element in eq.(|l]). If one could prove that such vacuum effects indeed 
reduce the value of G\ from 0.579 to 0.29, then, on the one hand, the 'proton spin crisis' 
would be resolved since the naive parton interpretation and the OZI estimate of G\ would 
be no longer valid. On the other hand, such a result would establish the existence of 
nontrivial topological vacuum structures, thus opening the door for a quantitative study 



of the connection between vacuum topology and instantons|10 . 

The lattice method provides an ideal tool to calculate the matrix element of eq.(^) non- 
perturbatively. Within this method two complementary approaches have been developed 
over the recent years. The first is to evaluate the connected and disconnected contributions 
to the flavor singlet matrix element directly, which are schematically depicted in fig.|l|. This 
"direct method" will be dealt with in the present paper. 

Alternatively, one can exploit the anomalous divergence equation of the flavor singlet 
axial vector current. Using this approach, the singlet coupling G\ is extracted from the 



correlation function of the proton propagator and the topological charge density]!]], [T2], |T5 |. 
Although this "topological method" circumvents the computationally expensive calculation 
of disconnected (fermionic) contributions, it has some intrinsic difficulties. The connection 
between G\ and the correlation of proton and topological charge density is valid in the 
limit of vanishing quark mass and zero momentum transfer only. As both these limits 
cannot be taken by a direct calculation on a finite lattice, two (wide range) extrapolations 
have to be performed. In addition, the "topological method" is ill defined in the quenched 
approximation of QCD|]TT|. On the other hand, the latter method provides the advantage 



of a nonperturbative renormalization procedure for the topological charge density |14, |15 



The "direct method" has been tested on quenched configurations by Fukugita et al. ]TB 
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and Dong et al. [|HJ some time ago. Although the former study has been done at a somewhat 
strong coupling {(3 = 5.7) and the latter is based on a very limited statistical sample (24 
gauge configurations), the results are encouraging, as they find estimates for G\ compatible 
with eq.(|]). 

It is of course by no means obvious, that these estimates hold also in full QCD. As we 
explained above, one expects vacuum contributions to lower the parton estimate of G\. 
The actual numbers might well depend on the details of the vacuum structure, such as the 
presence of sea quarks. Therefore, full QCD lattice simulations are necessary to calculate 
reliable QCD results for G\. 

Exploratory full QCD studies using the "topological method" have been performed by 
the authors of ref.|L3j and ref.[@, within the Kogut-Susskind lattice discretization scheme 
and with rif = 4 degenerate quark flavors. Both simulations suffer however from serious 
systematic uncertainties, as they have been performed at finite quark mass and sizeable 
momentum transfer (q 2 ~ (600MeV) 2 ). 

In this paper we present the results of a full QCD calculation from the "direct method" . 

We have analyzed 200 gauge configurations at each of our four sea quark masses cor- 
responding to m^/rrip = 0.833(3), 0.809(15), 0.758(11) and 0.686(11). The gauge configu- 
rations have been generated previously in the standard Wilson discretization scheme with 
n/ = 2 mass degenerate quark flavors, at (5 = 5.6 and with a lattice size of n^xrit = 16 3 x 32 
points. Details of the simulation can be found in ref. 1^ . 

As a preparing remark we emphasize that the numerical calculation of the disconnected 
parts of axial vector quantities is difficult, especially in full QCD simulations. It is known 
from lattice evaluations of the pion-nucleon-sigma term [^0], that the signals for the 
disconnected contributions of scalar insertions suffer from large statistical fluctuations. 
Since the disconnected parts of axial vector quantities can be seen as the difference of two 
such scalar insertions, one can expect that the statistical noise is even more dangerous in 
the axial vector case. On top of this, a comparison between quenched and full QCD lattice 
results for the light spectrum revealed [[HJ, that the statistical uncertainty in full QCD is 
larger by about a factor of 2 compared to a quenched calculation at equal statistics, cutoff 
and lattice size. Thus, one major issue of this paper is to investigate whether and to what 
accuracy one can extract disconnected signals for axial vector quantities with state of the 
art methods and statistics in full QCD calculations. 

The paper is organized as follows. In the next section we will briefly explain the methods 
used to calculate the connected and disconnected contributions to the matrix element of 
eq.(|l]) and present our raw data. The chiral extrapolations and the renormalization are 
performed in section 3. A discussion of the physics results and concluding remarks are 
given in section 4. 



2 Lattice Methods and Raw Data 
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2.1 Connected contributions 



To calculate the connected contributions to the axial vector density matrix element of the 
proton, eq. ([!]), we have applied the global summation method |2l], ^2fl . With this technique 
one calculates the ratio 



jdSUM 



(t) = 

E^ f (0, 0) Ey, yo frvYsg] % Vo)P(x, t)) 

E^(Pt(o,o)P(f,t)) 



(3) 



£ [flvrsg] (y,yo)) 



where P is an interpolating operator for the proton. From this ratio one extracts the 
matrix element (-P| 97^75? |-P), which is given by the slope of the asymptotic form , 



jdSUM 



W A+{N\q 1 ^q\N)t 



(4) 



Note that there are two different types of connected insertions, one from the contraction 
of the "loop" quark q with the it-type quark of the proton and one from the contraction 
with the d type quark. This yields 

Rq,u(t) = A q + C q ^t , C q ^ = (P(K sea ) \ [g7M75<?] sea)\P sea)) con , q = U,d . (5) 

The 3-point correlator in the numerator of eq.(||) has been evaluated using the standard 
insertion technique ||22||. To improve on the ground state projection of the proton we have 
applied 'Wuppertal smearing' |23[ with n = 50, a = 4 to all proton operators. 

In fig.|| we display the corresponding signals, spin averaged over 7^75, /i = 1, 2, 3. The 
sea quark mass has been set equal to the (proton) valence quark mass. We emphasize that 
the expected linear rise in i?^ c/M (t) is manifest for all four quark masses. 

The results of the fits to R^ M according to eq. (|5|) are collected in Table [I]. We 
find only a weak dependence on the quark mass. Note that the statistical errors to the 







C u 


c d 


0.1560 
0.1565 
0.1570 
0.1575 


0.833(5) 
0.809(15) 
0.758(11) 
0.686(11) 


1.100(8) 
1.102(11) 
1.025(12) 
1.018(24) 


-0.307(3) 
-0.308(4) 
-0.297(8) 
-0.295(10) 



Table 1: Lattice results for the connected spin averaged amplitudes C u and Cd at finite 
quark mass. 
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Figure 2: Summation method: The raw data R u and R,i for the spin averaged connected 
contributions at our sea quark masses. The fits (range and value) are indicated by dashed 
lines. 



2.2 Disconnected contributions 

The disconnected contributions to the axial vector density matrix element can be analyzed 
in principle with the global summation method too. It has been pointed out however in 
ref.pl] that the application of this technique induces large statistical fluctuations. We 
have therefore chosen a slightly modified procedure, the plateau accumulation method 
(PAM) I^TJ , which leads to a much better signal to noise ratio. The PAM ratio is defined 
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as 



t-At 

R p A AM (t,At ,At) = £ RS^&yo), (6) 

2/o=At 



where 



W>) = (7) 

£*(-P + (o,o)£7[ffVM] (y,yo)P(x,t)) „ 

(El^TMl to-*))- 

At and At can be varied in the range 1 < At, At < t to study the influence of contribu- 
tions from excited proton states. 

The asymptotic time dependence of R PAM is given by 

R P 4 M {t, At , At) = B + D^(t - At - At ) , (8) 

where D q denotes the disconnected contribution to the axial vector amplitude 

D q4l = (N\q W q\N) disc . (9) 

The determination of the numerator of eq.(^) requires the calculation of the correlation 
of the proton propagator and the axial vector quark loop La Avo) — Tr{q r )^ r ) 5 q) at a 
given timeslice y . La{vo) has been estimated using the spin explicit stochastic estimator 
technique |21| with complex Z2 noise |24| and 100 stochastic estimates per spin component 



and configuration. 

We display in fig.|]the signals[| R PAM (t), withQ At = At = 1, for the symmetric case, 
where the (proton) valence quark mass and the "loop" quark mass have been set equal to 
the sea quark mass. Note that the statistical errors are quite large. In particular for our 
lightest quark mass, corresponding to n sea = 0.1575, the signal is seen to be weak. We 
emphasize however, that the slope of R PAM appears not to depend strongly on the quark 
mass. Thus, excluding the result for the lightest quark mass from the analysis would not 
change the extrapolation to the chiral limit within errors. 

A similar behavior is found for the nonsymmetric case, where we have set the 'loop' 
quark mass equal to the strange quark massjl9|, while keeping the valence quark mass at 



the sea quark mass. We show the corresponding signals in fig. |] for the heaviest and the 
lightest quark mass. Within sizeable statistical errors one can identify the expected linear 
decrease reasonably wellQ. 

In Table ^| we have collected the raw results for D q 3 and D s 3 at the different quark 
masses. Apparently, the disconnected contribution to G\ is definitely less than zero, albeit 
within statistical errors of the order of 50%. 



1 Note that we have used only the combination 7375 in eq. (0), since we found the signals corresponding 
to the combinations 7175 and 7275 to be dominated by statistical noise. 

2 The signals R^ AM (t) with Ato = At > 1 are too noisy to extract a reliable signal. 

3 The same situation is found at K sea — 0.1565. For K sea = 0.1570 however, we observed a linear 
decrease for small t, though the asymptotic behavior is not clear. 
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Figure 3: The ratio R^ M (t) with Ato = At = 1 for the disconnected symmetric ampli- 
tudes D qt3 = (P(Ksea)\{qj3l5<l)(K sea )\P(K sea )) at our sea quark masses. The fits (range and 
value) are indicated by solid lines. 



3 Results 

3.1 Chiral Extrapolations 

In order to arrive at our lattice estimate for G\, we have to extrapolate the connected 
(C u , Cd) and the disconnected amplitudes (D q , D s ) to the light quark mass. Since the 
"loop" quark mass in D s is kept fixed at the strange quark mass, this quantity serves as a 
semi quenched estimate for the strange quark contribution to G\. 

We display the (unrenormalized) lattice data and the corresponding fits as a function 
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Figure 4: The ratio R^? M {t) with Ato = At = 1 for the disconnected non symmetric 
amplitudes D St3 = (P(K se a)| (973759) (fts)l-P(ftsea)) at the heaviest and the lightest sea quark 
mass, and k s = 0.15608. The fits (range and value) are indicated by solid lines. 



n 


m n /m p 


D q 


D s 


0.1560 
0.1565 
0.1570 
0.1575 


0.833(5) 
0.809(15) 
0.758(11) 
0.686(11) 


-0.146(60) 
-0.177(97) 
-0.119(62) 
-0.130(69) 


-0.163(58) 
-0.125(74) 
-0.129(82) 
-0.184(76) 



Table 2: Lattice results for the disconnected amplitudes D q3 = (^(k)] [973759] (ft)|-P(fv))<2i.; 
and D s 3 = (P(n)\[qj 3 j 5 q](K s )\P(K,)) ( ns- The value of the hopping parameter k s = 0.15608 



corresponds to the strange quark mass ]7^/ 



of the (sea) quark mass in fig.|5|. As mentioned above, the mass dependence of all the 
amplitudes is weak and the linear ansatz for the extrapolations is in accord with our data. 
Within errors, D q and D s do agree. We conclude that, to the present level of accuracy, 
flavor symmetry of the disconnected parts appears to be maintained in (nf = 2) full QCD. 
Note that the disconnected contribution to G\, given by 2D q + D s , is larger than the 
connected contribution Cj, suggesting a sizeable modification to the OZI inspired estimate 
of the singlet coupling. 

At the light quark mass, being determined by the requirement that the ratio m^/irip 
on the lattice equals the experimental result we find 

C u = 0.953(25) , C d = -0.290(11) , 
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Figure 5: Extrapolation of the unrenormalized lattice amplitudes to the light quark mass. 
The fits are indicated by dashed lines, the results of the fits by bursts. To avoid cluttering 
the data D s (m q ) have been slightly shifted. 



(10) 

D q = -0.108(63) , D s = -0.166(91) . 

The numbers in brackets denote the statistical errors, obtained by a jackknife analysis of 
the extrapolation. 

From the appropriate combinations of these amplitudes one can extract the lattice 
values of the singlet and nonsinglet couplings, as well as the contributions of each single 
quark flavor to these quantities. The results are collected in Table |3|. Note that the 
nonsinglet couplings have been calculated assuming flavor symmetry of the disconnected 
parts. 
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Au = C u + D q 


A d = C d + D q 


A s = D s 


0.85(10) 


-0.398(97) 


-0.166(91) 


G\/Z S A = C U + C d + 2D q + D s 


ri'A IvNS _ r 1 r> 
(j ' Al ^ 'A ~ w — W 


r-<8 17NS — n j_r> 
^ Al ^ A —^u-r^d 


0.28(16) 


1.243(28) 


0.663(25) 



Table 3: Lattice results at the light quark mass for the contribution of each quark flavor 
and for the axial vector couplings of the nucleon. (Z A S ),Z A denotes the (non)singlet axial 
vector renormalization constant . 



3.2 Renormalization and Physics Results 

In continuum QCD the flavor nonsinglet current A^ s is conserved in the chiral limit. Thus 
it needs no renormalization. The flavor singlet axial vector current on the other hand 
is not conserved due to the presence of the axial anomaly. This leads in the continuum 
to a nontrivial renormalization of G\. In second order perturbation theory picks up 
a logarithmic divergence in the cutoff, which induces a renormalization scale dependence 
on the axial flavor singlet coupling, G\ = G\(/i 2 ). The associated anomalous dimension 
is known to 3 loops|2^]. Fortunately, it turns out that the /i 2 dependence of G\ is very 
weak, e.g. increasing fx 2 from lOGeV 2 to infinity decreases the value of G\ by about 10% 
only, see ref. ||. 



The situation is slightly more involved in the Wilson discretization of lattice QCD |26| . 
Already the nonsinglet current A^ s is not conserved and has to be renormalized by a finite 
factor Z^ s . On top of this the singlet current suffers from the axial anomaly. The resulting 
scale dependence of Z\ might be different from the continuum form, as the axial anomaly 
mixes with the Wilson term. 

In full QCD both, Z^ s and Z A , are known only to first order (tadpole improved) 
lattice perturbation theory |27|, [28[] , where they agree. For the nonsinglet current, non- 
perturbative renormalization procedures have been developed and used in the quenched 



approximation 1 29, 30]. Those methods, when applied to the case of full QCD, will yield a 



NS 
A 



reliable value for Z 

The /1 2 dependence of Z A however, which occurs beyond first order, has not been 
calculated yet. One may hope of course, that it is as weak as in the continuum. 

Thus, in view of the current state of analysis, we have used the first order tadpole 

7S 

J A-. 



improved estimate for Z A = Z'J 



with ajfsd) = 0.215, k c = 0.158507, and k = «, = 0.15846, k = k s = 0.15608 for the light 
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and strange quark insertions respectively |]T9 . 

Unfortunately this choice makes it difficult to determine the scale /z 2 at which we 
calculate G\. The renormalized results are compiled in Table |j. Comparing the estimates 



Au 


A, 


A s 


G A 


G\ 




0.62(7) 


-0.29(6) 


-0.12(7) 


0.20(12) 


0.907(20) 


0.484(18) 



Table 4: Renormalized results at the light quark mass for the contribution of each quark 
flavor and for the axial vector couplings of the nucleon. For the definitions of the quantities 
see Table | . 



for G\ and G\ one finds that the disconnected contributions lead to a violation of the 
Ellis- Jaffe sum rule by more than 50%, which is consistent with the experimental finding. 
Note that there is a 30% discrepancy between our estimate of the triplet coupling G A and 
the experimental value, G\ = 1.2670 ± 0.0035|3l|. This points to the presence of sizeable 



higher order or even nonperturbative contributions to Z A S . An uncertainty of similar size 
in Z\ would however be well covered by the statistical errors of 50% in G\. 



4 Discussion and Outlook 

We have calculated connected and disconnected contributions to the flavor singlet axial 
vector coupling of the proton in a full QCD nf = 2 lattice simulation with Wilson fermions. 
We find G\, within large uncertainties, to be consistent with the experimental result and 
with previous quenched estimates. This indicates that the substantial violation of the 
Ellis- Jaffe sum rule might be caused mainly by the gluonic properties of the vacuum. 

The major sources of uncertainty are the large statistical fluctuations of the discon- 
nected parts and the use of the 1st order estimate for Z A . 

The statistical quality of disconnected axial vector insertions, achieved with state of 
the art stochastic estimator methods and with a sample size of 200 configurations, is 
encouraging but not yet satisfactory. Ideas to improve on this situation are to use larger 
lattices, where self averaging effects will help, to adapt the stochastic estimator techniques 
to the axial vector case, e.g. by use of correlated noise and, last not least, by an increased 
statistics. Work along these lines is in progress. 

The nonperturbative determination of Z A presents a major problem. One possible 
solution is to perform a nonperturbative renormalization of the topological charge and then 
to use the 'topological method' to extract the flavor singlet axial vector renormalization 
constant. 

Having overcome these problems one has to study cutoff and volume dependence of G A . 
Finally, as a future project, one should perform full QCD lattice simulations with nj > 3 
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non degenerate dynamical fermions. This would allow for a realistic estimate of D s and of 
the amount of flavor symmetry breaking of the disconnected contributions. 
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